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To investigate the non-perturbative, electric sector of a deconfined gauge theory at nonzero tem¬ 
perature, we consider a SU( 2) matrix model. We compute /3-functions to one loop order for the 
simplest extension of the 0(4) nonlinear sigma model, which involves three coupling constants. 
Computing in the ultraviolet limit in 2 + e dimensions, we find that at least one coupling is not 
asymptotically free. 

PACS numbers: 


I. INTRODUCTION 

The deconfining phase transition in SU(N) gauge the¬ 
ories is of interest in its own right, as a problem in sta¬ 
tistical mechanics, and for its possible application to the 
collisions of large nuclei at very high energies. 

If the deconfining transition is of second order, then 
universality and the renormalization group predicts the 
critical behavior [1], Lattice simulations indicate that 
the deconfining transition is of second order for TV = 2, 
N = 3, and (perhaps) N = 4 in 2 + 1 dimensions [2], and 
for TV = 2 in 3 + 1 dimensions [3]. 

For other N, N > 5 in 2+1 dimensions [2], and for TV > 
3 in 3+1 dimensions [4, 5, 6, 7], the deconfining transition 
appears to be of first order [4, 5, 6, 7]. Universality is of 
no help for a first order transition. Moreover, even when 
a transition is of second order, universality is only of use 
near the transition; it says nothing about what happens 
away from it. 

A first order transition is natural at large N , because 
the number of degrees of freedom is ~ 1 in the confined 
phase, and ~ N 2 in the deconfined phase [7]. One pos¬ 
sibility is that the transition is so strongly first order 
that the entire deconfined phase is a nearly ideal gluon 
plasma. In 3+1 dimensions, for example, this might have 
been true for all TV, from TV = 3 to N = oo. In this case, 
the deconfined phase would be amenable to analysis by 
means of a resummed perturbation theory [8], for all tem¬ 
peratures from the transition temperature, on up. This 
does not seem to be what happens, however. If Ty is the 
temperature for deconfinement, all consistent resumma¬ 
tions of perturbation theory appear to fail in describing 
temperatures below temperatures of order ~ 4Ty [8] . 

A plausible guess is that at temperatures T between 
Td and ~ 4 Td, the theory goes into a regime of strong 
coupling [9, 10]. After all, by running asymptotic free¬ 
dom in reverse, as the temperature decreases, the value 
of the strong coupling constant, a s {T ), increases. 

Even at temperatures as low as the transition temper¬ 
ature, though, the gauge coupling is not especially large. 
Consider the effective coupling which enters into the 


dimensionally reduced theory in three dimensions [11]. 
Computations show that in the magnetic sector, noth¬ 
ing surprising happens, even at temperatures as low as 
Td- For example, one can compare the static string ten¬ 
sion in the dimensionally reduced theory, to that in the 
full theory. Using the two loop calculations of the mag¬ 
netic coupling by Giovannangeli [12], Laine and Schroder 
[13] find remarkably good agreement between these two 
quantities from high temperature all of the way down to 
Td- They estimate that in QCD, even at a “transition” 
temperature of « 175 MeV, that the effective coupling 
in the dimensionally reduced theory is ss 0.28 [13]. 
This is a surprisingly small coupling: at zero tempera¬ 
ture, it corresponds to a relatively high momentum scale, 
of « 1.6 GeV. 

While nothing remarkable happens in the magnetic 
sector, something striking happens in the electric. At 
nonzero temperature, the electric sector of a gauge the¬ 
ory is probed by the eigenvalues of the thermal Wilson 
line, which are gauge invariant [14, 15, 16, 17, 18, 19, 20, 
21, 22, 23]. Most notably, the Polyakov line is propor¬ 
tional to the sum of these eigenvalues, and is near one 
in a perturbative regime. Numerical simulations on the 
lattice show that while the renormalized Polyakov loop 
is ~ 0.9 at « 3 Td, it falls sharply as the temperature 
decreases, to ~ 0.4 at the transition, and vanishes in the 
confined phase [19, 20]. 

Thus when a gauge theory deconfines, it does not go 
immediately to a nearly ideal plasma at Td- Instead, for 
temperatures between Td , and « 4 Td, the electric sector, 
although not strongly coupled, behaves in a pronounced 
non-perturbative manner. 

We remark that this behavior in the electric sector 
drives the theory far from the conformal limit. This is 
shown by the interaction measure, which is the trace of 
the energy momentum tensor divided by T 4 . The in¬ 
teraction measure vanishes for a conformally invariant 
theory, such as for a free, massless field, or less trivially, 
for a gauge theory with Af = 4 supersymmetry. For an 
SU(N) gauge theory in 3 + 1 dimensions, lattice sim¬ 
ulations show that for TV = 2 [3], TV = 3 [4], TV = 4 
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[6], and TV = 8 [6] , the interaction measure is very large 
when T : Td —> 4T^. How large can be estimated by 
comparing to the value in perturbation theory, where it 
starts at ~ For 77 = 3, the interaction measure has a 
sharp maximum just above the transition, at T « 1.2 Tj,- 
At this maximum, the interaction measure is about ten 
times larger than its perturbative value, using the values 
of u s (T) from [13]. In contrast, by « 4 Tj, its value is 
within the range expected from perturbation theory. 

In order to understand the non-perturbative, electric 
sector of a deconfined gauge theory, it is necessary to 
develop an effective theory of thermal Wilson lines [16, 
20, 21, 22, 23]. Wilson lines are SU(N) matrices, and so 
the effective theories of relevance are matrix models. 

In this paper we take a small step in this process, by 
considering the effective theory for deconfinement in 2 + 1 
dimensions. In 3 + 1 dimensions, at least at large N 
the effective theory for deconfinement is dominated by 
an effective potential [22]. In 2 + 1 dimensions, however, 
fluctuations can dominate the potential. Fluctuations are 
why, for instance, the deconfining transition for N = 3 is 
of second order in 2 +1 dimensions, versus the first order 
transition expected in mean field theory. 

The infrared behavior of such a matrix model is in¬ 
volved, and involves complex interactions between the 
potential and fluctuations. The ultraviolet limit, though, 
is amenable to perturbative analysis, and it is this which 
we study in the present work. 

Our computation is similar to that for an ordinary non¬ 
linear sigma model in two dimensions [24, 25, 26, 27, 28, 
29] . While familiar sigma models have one coupling con¬ 
stant, matrix models have an infinity of couplings, all 
of which are relevant in the ultraviolet limit. Thus one 
of the purposes of this paper is to see if such models 
are well behaved in perturbation theory. We find that 
they are: modulo a reasonable technical assumption on 
the size of the coupling constants, (48), we find that at 
least to one loop order, the counterterms induced can 
be reabsorbed into renormalizations of the bare coupling 
constants. This allows /9-functions to be defined in a 
standard fashion [24]. 

It is well known that the single coupling constant of 
the nonlinear sigma model is asymptotically free in two 
dimensions [24, 25, 26, 27, 28, 29]. One might hope that 
the new couplings in matrix models are also asymptot¬ 
ically free. We find that this is not true. We consider 
the simplest generalization of an 517(2) x 517(2) sigma 
model, which involves three coupling constants. While 
the coupling analogous to that of the usual sigma model 
is always asymptotically free, at least one of the two new 
couplings is not. This result is similar to that of Friedan, 
who considered a sigma model in a background metric, 
and found that the new couplings induced by the back¬ 
ground metric usually spoil asymptotic freedom [27]. 

In Sec. II we classify the types of matrix models for 
general N; in Sec. Ill, for N = 2. Sec. IV describes 
the one loop calculation in a general background field. 
The counterterms are computed in Sec. V, producing 


the /3-functions of Sec. VI. In Sec. VII we suggest how 
the /3-functions for N = 2 might generalize to 77 > 3, 
and how these affect the phase transitions of the model 
in, and above, two dimensions. An appendix contains a 
comment about a nearly Goldstone boson at large N in 
SU(N) gauge theories. 

II. VARIETY OF MATRIX MODELS 

Consider a matrix valued field L(£). Let it be a SU ( N ) 
matrix, satisfying the constraints 

L^L = 1 , det(L) = 1 . (1) 

We chose L to lie in the fundamental representation, so 
L = L ab, where a, b = 1... 77. In the usual nonlinear 
sigma model, this matrix is assumed to be invariant un¬ 
der a global symmetry of ( SU(N ) x SU(N))/Z(N), 

L{x) -» U L(£) V ; (2) 

U and V are independent global 517(77) rotations, mod¬ 
ulo a common Z(N) rotation. If U and V are distinct 
rotations, the only renormalizable Lagrangian invariant 
under this symmetry is 

-4 tr l 5 * L | 2 > ( 3 ) 
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where g 2 is the coupling constant for the sigma model. 

For a generic sigma model in two spacetime dimen¬ 
sions, both the field, L, and the coupling constant, g 2 , 
are dimensionless. It is well known that the coupling g 2 
is asymptotically free in two space-time dimensions [24]. 

The model can be generalized by relaxing the global 
symmetry to one of SU(N) J taking V = W. The com¬ 
mon assumption is to impose a further constraint on the 
trace of L. Then the only possible action remains as 
in (3), but the symmetry changes, as the model now 
defines a symmetric space [24, 25, 29]. For example, 
if 77 = 2 M is even, and one imposes the constraint 
that L is traceless, tr L = 0, the symmetry becomes 
SU(2M)/S(U(M) x U(M)) [24, 25, 29], 

To study the deconfining phase transition in d + 1 
spacetime dimensions, one can construct an effective the¬ 
ory of straight, thermal Wilson lines in d spatial dimen¬ 
sions [23]. This effective theory is valid over (spatial) 
distances x /§> 1/T, where T is the temperature. As 
a theory in two spatial dimensions, then, the model we 
study is relevant to the transition in 2 + 1 dimensions. 
We will discuss the salient properties of this effective the¬ 
ory as we go along. For now, we only need to know that 
in the effective theory, the Wilson line, which we denote 
as L(x), is invariant under local gauge transformations, 
f2(x), and global Z(N) transformations, 

L(£) -+ e 2 ™/ N H(£) L(£) fV(£) . (4) 

To accommodate the local gauge invariance, it is nec¬ 
essary to include a SU(N) gauge field, A t . Define the 
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adjoint covariant derivative as A = di + iG[Ai,], where 
G is the gauge coupling. The action for L is then 

tr |AL| 2 , (5) 

plus the usual action for the gauge field. In the effective 
theory for deconfinement, the A,j represent the gauge po¬ 
tentials for the (static) magnetic field. 

Unlike symmetric spaces, the trace of L cannot be con¬ 
strained. In fact, the simplest trace, 

£ = h tr L ’ ( 6 ) 

is an order parameter for deconfinement, the Polyakov 
loop in the fundamental representation. A nonzero value 
for the fundamental loop signals that the vacuum spon¬ 
taneously breaks the global Z(N) symmetry in the de- 
confined phase. 

As we cannot impose a constraint on the trace of L, 
the possible Lagrangians are far more complicated than 
those of the usual nonlinear sigma model. To start with, 
there are terms with no derivatives; that is, there is a 
potential for L. This is a sum over loops [22], 

V= X] «5 Re Zs ■ (7) 

SelZo 

Here, is denotes the trace of a loop in the representation 
S. For example, the loop in the adjoint representation is 

iad = ( l tr L ! 2 “ X ) ■ ( 8 ) 

We always divide a loop by the dimensionality of the 
representation. In this way, in the perturbative limit, 
when L = 1, all loops are normalized to be one. 

We assume that the breaking of the global 
Z(N) symmetry is spontaneous, so that the only 
loops which contribute to the loop potential are 
Z{N) neutral. In the loop potential, TZq denotes all possi¬ 
ble Z{N) neutral representations. The series starts with 
the adjoint loop. Using the character expansion, we only 
need take linear powers of loops, albeit in arbitrary rep¬ 
resentations. 

Ref. [23] discusses the construction of an effective the¬ 
ory in 3 + 1 dimensions, but at least formally, it is im¬ 
mediate to extend it to 2 + 1 dimensions. In any number 
of dimensions, the Wilson line is dimensionless. For the 
original gauge theory in 2 + 1 dimensions, the gauge cou¬ 
pling squared, g 2 , has dimensions of mass. In the effective 
theory, classically the term for the electric field is identi¬ 
cal to (3), with the sigma model coupling, g 2 = g 2 /T [23]. 
The gauge coupling in the effective theory is G 2 = g 2 T, 
and so G 2 has dimensions of mass squared. At one loop 
order, some potential terms, as in (7), are induced, with 
ks ~ T 2 . Further terms are induced by corrections to 
higher order; these are then a power series in g^/T, times 
T 2 . Besides terms with no and two derivatives, terms 


with four and more derivatives also arise in the effective 
theory [23] . We stress that the construction of this effec¬ 
tive theory is only valid in perturbation theory, where 
g 2 /T is small; implicitly, it is assumed that this is a 
reasonable approximation, even down to the transition 
temperature. While this appears to be true in 3 + 1 di¬ 
mensions [13], this has not yet been studied in lattice 
simulations in 2 + 1 dimensions [2] . 

The dominant behavior of the theory in the infrared 
and ultraviolet limits can be read off by the dimension¬ 
ality of the couplings. In the infrared limit, as all ks 
have dimensions of mass, they should dominate. There 
is an important, if familiar qualification: as a theory in 
two dimensions, fluctuations can be important when the 
fields are light. 

Conversely, in the ultraviolet limit, terms with two 
derivatives dominate over those with none. Now of 
course, since we are ultimately interested in an effective 
theory, valid only over large distances, studying the ul¬ 
traviolet behavior is something of an academic exercise. 
It was our original hope that these theories might be 
asymptotically free in all couplings. If so, then as with 
the usual nonlinear sigma model, one could be able to 
study the phase transition in the infrared limmit, in three 
dimensions, by developing an expansion in 2 + e dimen¬ 
sions, e > 0. In the next section and the remainder of 
the paper, we consider the simplest extension of the non¬ 
linear sigma model, which involves the three couplings 
which contribute to quartic interactions in the perturba¬ 
tive limit, Sec. IIIB. While that leaves an infinity of 
other couplings, since we find that these three couplings 
are not uniformly asymptotically free, Sec. VI, there is 
no point in going any further. 

With this mea culpa aside, we classify all terms with 
two derivatives [22], The most obvious is the original 
Lagrangian, times any Z(N) neutral loop: 

tr |AL| 2 Re A , es = 0 . (9) 

Re denotes the real part, while es is the Z(N) charge of 
a loop in the representation S. As Z(N) charges, es is 
only defined modulo N. For the time being, we do not 
bother with denoting coupling constants, since we are 
only concerned with the types of terms which can arise 
in the action. 

One can add more powers of L inside the trace in (9), 
Retr (|AL| 2 L p ) i s , e s = -p . (10) 

The global symmetry is only one of Z[N) , and not 
C7(l). This implies that instead of taking the trace of the 
complex conjugate square of the covariant derivative, as 
in (9), one can take the trace just of the square of the 
covariant derivative: 

Re tr (AL) 2 l S , e s = -2 ; (11) 

this is then multiplied by a loop of charge —2 to ensure 
Z(N) invariance. Further powers of L can be added in- 
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side the trace: 

Retr ((AL) 2 I/) £ s , e s = -p - 2 . (12) 

This list continues, until one is only left with deriva¬ 
tives of loops: 

Re (diis ) (di£s>) £s n , es + es> + es n = 0 . (13) 

This list is not meant to be exhaustive. Clearly there 
is an infinite set of relevant couplings. As we shall see 
in the next section for N = 2, not all of these couplings 
are necessarily independent. While this infinite set of 
couplings is much more complicated than ordinary sigma 
models [24, 25, 26, 27, 27, 28, 29], it is dictated by the 
physics, and is unavoidable. 

For models with SU(N) x SU(N) symmetry, and their 
associated symmetric spaces, the exact 5-matrix can be 
determined by Bethe ansatz [28]. Because of the plethora 
of couplings, we do not expect that our matrix model is 
soluble by similar means. 

III. N=2 MATRIX MODELS 
A. Classifying Lagrangians 

We turn to the case of two colors. Two colors is illu¬ 
minating, since one can explicitly compute all elements 
of the Lie group using the exponential parametrization: 

L = exp (i a -t) , (14) 



We use generators t a , normalizing tr (t a t b ) = S ab /2 (so 
t a = cr a / 2, where the o a are the Pauli matrices). We also 
denote 

a = a a , a 2 = 1 . (15) 

The Lagrangian of the gauged nonlinear sigma model 
is 

£o(L) = 2 tr | AL| 2 . (16) 

We defer the definition of coupling constants until later. 
For N = 2, the covariant derivative is 

A a = di a + G A t x a . (17) 

The gauge field only couples to the isovector, a a , and 
not to the isoscalar, a. For N = 2, the Lagrangian of the 
sigma model is then 

C 0 {a) = (di a) 2 + 4 sin 2 (AS f . (18) 

To avoid clutter, the dependence of £o(a) on the gauge 
field, is left implicit. 


The next term, as in the series in (10), is to multiply 
the adjoint loop times this term. We subtract off the 


value of the loop when L = 1, and so define 

A(L) = | (l-4d) tr |AL| 2 , (19) 

= sin 2 £ 0 (S) . (20) 

Another possible kinetic term involves the loop in the 
fundamental, or doublet, representation: 

£ 2 (L) = 4(c^) 2 , (21) 

= sin 2 (did f . (22) 


This is a term as in (13). As this is formed exclusively 
from the loop, it depends only upon the magnitude of 
d?, and not upon its direction in isospin space, a. For 
SU( 2), £ = tr L/2 is automatically real, so it doesn’t 
matter if one takes just the square of the derivative of £, 
or the complex conjugate square. Lastly, as the loop is 
gauge invariant, A is independent of the gauge field. 

A term like that in (12) is 

tr (AL) 2 = -Co(S) + A(o0 - (23) 

This identity is special to N = 2, and shows that this is 
not a new, independent coupling. 

It is clear that there is an infinity of possible cou¬ 
plings. In going from the original Lagrangian, (16), to 
(19), we multiplied by the adjoint loop minus one, which 
for TV = 2 is proportional to sin 2 (a/2). We can continue 
this process to infinite order, multiplying A and £2 by 
higher and higher powers of sin 2 (a/2). All of these are 
independent couplings, with dimensionless coupling con¬ 
stants. 

For N = 2, the Lagrangian C 0 is invariant under an ex¬ 
tended global symmetry of (SU( 2) x SU(2))/Z(2), which 
is isomorphic to 0(4) [26]. It is instructive to write the 
Lagrangians in terms of 0(4) fields. This is done by in¬ 
troducing 

er = cos(a/2), 7? = sin(a/2)d , (24) 

and then forming the 0(4) vector </> = (a, if). The matrix 
constraint then becomes </ 2 = ct 2 +7? 2 = 1. Dropping 
the gauge fields, we find 

C 0 ($) = 4 (di$) 2 = 4 ^(d ; ;7r) 2 + ^ ^ . (25) 

The doublet loop i = a, so the new couplings in our 
matrix model are 

£iA) = (1 - a 2 )(di$) 2 = if 2 ((diH) 2 + ^ ^ , 

(26) 

and 

£ 2 (</)=4(cA) 2 =4 ■ (27) 

1 — 7 T z 
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It is clear that only Cq is 0(4) invariant. We write C\ 
and £2 using both 0(4) degrees of freedom, and in terms 
of the 7? field. We do this because while the independent 
degrees of freedom are the tt fields, it is instructive to see 
how the 0(4) symmetry is broken in the new couplings. 
In terms of the ff’s, (23) is especially simple, 

tr (A:L) 2 = - 4 (d^) 2 , (28) 

just a free kinetic term for the pions. 

The /3-functions can be computed using the 7r’s, but 
we find the a’s more convenient. 


B. Coupling Constants 

We introduce the coupling constants g 2 , ft and A in 
the bare Lagrangian as 

£ = ~z 2 £o + 4 £i + — £2 ■ (29) 

2 g 2 g 4 g 4 

To understand why we introduce the coupling constants, 
£ and A, with overall factors of 1/g 4 , we first note that 

£°(a) = ^ sin 2 (^) ((Did) 2 - (dio) 2 ) + (dia) 2 . (30) 

This form is useful in the perturbative limit. The sim¬ 
plest perturbation theory is to expand about zero held, 
taking a => g/3. Expanding to quartic order in /3, 

±-Co = \(D i pj 2 -^ 1 p 2 (DS) 2 +... , (31) 

4£i « j/3 2 (A/3) 2 +... , (32) 

5 4 

and 

4 £2 « 4/3-ft/3) 2 +... . (33) 

9 4 

Thus £0 includes the free part of the Lagrangian, plus 
quartic interactions ~ g 2 . For zero background held, the 
new terms in the Lagrangian, C\ and £ 2 , do not con¬ 
tribute to the quadratic part of the Lagrangian: they 
only contribute to quartic interactions ~ £ and ~ A, 
respectively. Hence the couplings g 2 , ft and A are all 
couplings which should be included to one loop order. 

This continues with the terms neglected above: a term 
such as (sin(/3/2)) 2 "£ 0 hrst contributes to a 2n + 2 point 
function of the /3’s. Thus the proper normalization of 
such a term is a new coupling constant times 1 /g 2(jl+1 \ 
This explains why we concentrate on the above three 
couplings: they are uniquely the only couplings which 
contribute to quartic interactions in the perturbative 
limit. 


IV. BACKGROUND FIELD METHOD 
A. Possible Classical Fields 

We wish to compute the renormalized Lagrangian to 
one loop order. As is well known, to do so it is easiest 
to use the background held method [24]. We take some 
classical, background held a, and expand in a quantum 
held (3. To one loop order, it is only necessary to expand 
to quadratic order in the quantum fluctuations, (3. The 
price paid is that all dependence upon the background 
held must be kept. Thus to ease our labor, we want to 
choose the simplest possible background held we can. 

The background held cannot be too simple, however. 
To one loop order, various counterterms are generated. 
We need a background held which allows us to distinguish 
which counterterms contribute to the renormalization of 
which terms in the bare Lagrangian. 

Consequently, even in principle, it does not suffice to 
expand about a trivial background: while the quartic 
interactions in £2 differ, (33), those in £0 and £ 1 , (31) 
and (32), are the same. 

Thus we need to expand about a background which is 
not trivial. The next simplest possibility is to take a held 
which lies everywhere in the same direction in isospin 
space, so that fta = 0. For such a held, the Lagrangian 
becomes 

£=M 1+ ^ sin2 ©) (fta)2 ' (34) 

This background held allows us to separate contributions 
to £0 from those to £1 and £ 2 , but we can’t disentangle 
which terms contribute to £ 1 , and which to £ 2 . In terms 
of f3- functions, we could determine that for g 2 , and the 
sum of ^ + A, but not for £ and A by themselves. 

Another possibility is to expand for a hxed direction in 
isospin space, but in the presence of a background gauge 
held, using the covariant derivative, [A, 3], to separate 
the various terms. We did not do this, because in order to 
respect gauge invariance, it is necessary to include fluc¬ 
tuations in the gauge held. (This was checked by explicit 
calculation.) Including the fluctuations in the gauge held 
seems unduly complicated, since they don’t contribute 
to the ultraviolet limit in two dimensions; their presence 
would be merely as a bookkeeping device, to sort out the 
different terms in the renormalized Lagrangian. 

The hnal alternative is to expand about a completely 
general background held, whose isospin direction changes 
in space, fta ^ 0. For an arbitrary held a, it is clear that 
£o(a?), £i(a), and £2 (a) represent distinct Lagrangians, 
(18), (20), and (22). Thus we can certainly pick out the 
contributions to different terms in the Lagrangian which 
are generated at one loop order. 

Indeed, the necessity of distinguishing between differ¬ 
ent terms in the renormalized Lagrangian is why we limit 
our calculations to two colors. The three couplings which 
we consider for N = 2 exist for N > 3, given in (16), (19), 
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and (21). Taking the classical field as L c i = exp (ia a t a ) 1 
with t a the generators of SU (AT), it is natural to take the 
a a to lie in the Cartan subalgebra. Indeed, the simplest 
form is to take a a proportional to a single generator, re¬ 
lated to global Z[N ) transformations, 


t • (a x / 3) = e abc t a a b /3 c , where e abc is the completely 
antisymmetric tensor. Further, 

L qu d t L f g U ss -i 80- t + ^t-0x 80 ) + ... (41) 


a a = a t N , t N =( _ (7V °_ 1} ) ■ (35) 

This is precisely like the SU(2) ansatz where dia = 0, 
(34), since in SU(2) we can always chose the fixed direc¬ 
tion in isospin space to lie along t 3 . As for N = 2, when 
N > 3 the ansatz of (35) is not sufficient to distinguish 
between £ and A, with the sum of £o + £i + £2 like that 
of (34). Thus to determine the separate /3-functions for 
£ and A, it is necessary to expand in a field which lies in 
two distinct directions. Since there are N — 1 elements of 
the Cartan subalgebra, when N > 3 we can chose the a a 
to lie in two, different elements of the Cartan subalgebra. 
Thus extending our calculation from N = 2, to N > 3, is 
straightforward to do. It is not difficult to see, however, 
that this will be tedious; even the form of the classical 
field is messy. In Sec. VII, we comment on results for 
N > 3, using the ansatz of (35). We use this to make an 
obvious guess about the leading form of the /3-functions 
when N > 3. 


To quadratic order, the doublet loop is 

= coe(|)-|sta(|)8-| 

-TT“ s (l)^ + - < 42 > 

Expanding the entire Lagrangian to linear order in 
the quantum fluctuations, and integrating derivatives by 
parts so that none acts on /3, we find 

6 X £ = - ■ (- diX (43) 

+ A. sing) (a 2 cos g)) a 

+ | g 0 (a)sin(a) S-4dgsin 2 g) a))^ ■ 
We introduce the modified couplings 


B. Explicit Expansion 




(44) 


We expand 

L = L ci L qu = exp(i a-t) exp (ig (3 ■ t ) , (36) 

where a is the classical, background field, and g [3 the 
quantum field. The factor of g is introduced to simplify 
later results, as in Sec. IIIB. This ansatz is convenient 
in computing £q(L) and £i(L), since 

tr (dilML) = tr^L^Lc (37) 

+ 2 L qu 5 s Lt u l4, diLct + d t L J u diLqu ) . 

Instead of the ansatz in (36), we could have taken L = 
exp(*(a + <7/3) -t). While this appears simpler, the expan¬ 
sion of (37) is then much more complicated. 

The product 

I 'h 8iL c i =i Ai-t, (38) 


as these arise naturally in the computation to one loop 
order. 

The equation of motion is 5 l C = 0. That of the ordi¬ 
nary nonlinear sigma model is just the first term on the 
left hand side in (43), diAi = 0. The additional terms, 
~ £ and A, are new contributions to the equation of mo¬ 
tion in this matrix model. 

The computation of terms to quadratic order in /3 is 
straightforward. We integrate by parts freely, organizing 
terms so that as few derivatives as possible act on /3’s, and 
as many on the background field, a. Doing so, we find 
that we can organize terms according to the maximum 
number of derivatives which act on /3. 

Terms with up to two derivatives acting on /3 are 

= 7/ (l- + 2^sin 2 g)) 0iP) 2 
+ A (di (sin g) a •/?)) . (45) 


is like a gauge field, albeit one formed from a pure gauge 
transformation, generated by the field L c j. Explicitly, 

Ai = dia a + sin(a) fta + 2 sin 2 g) a x dia , (39) 

where we have used the identity 

\ —t % — » —* 

a-t /3-f=-a-/31 + -f'(ax/3); (40) 


In writing <5|£, the superscript denotes an expansion to 
quadratic order in /3; the subscript denotes the maximum 
number of derivatives acting on /3. 

Terms with one derivative of /3 are 

«5 2 £ = g( 1 + 2esin 2 g)) A -0x80) 

+ £ sin a (a ■ (3) A • dip . 


(46) 
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Lastly, terms with no derivatives with respect to (3 are 

5lC = A cos(|) (<9 2 cos(^-)) /3 2 (47) 

+ | A) (a) (/3 2 - sin 2 (p 2 + (a • d} 2 )) . 

Because we have normalized the quantum field to in¬ 
clude a factor of g , terms to quadratic order which are 
independent of any couplings arise from the Lagrangian 
of the usual nonlinear sigma model. These are the first 
terms on the left hand side in (45) and (46). The 
first term on the left in (45) is just the usual free ki¬ 
netic term, (di/3) 2 / 2. The first term on the left in (46), 
—Ai ■ ( (3 x di/3)/2, represents the interactions in the usual 
nonlinear sigma model. All other terms, which are pro¬ 
portional to either £ or A, arise from the new couplings 
in a matrix model. 

There is an important restriction on the coupling con¬ 
stants in these models. Of course we can only compute 
in weak coupling, when g 2 , A, and £ are all <C 1. Since 
the new couplings in a matrix model affect the kinetic 
term, though, so that the usual kinetic term dominates 
at high momentum, it is necessary that the new terms 
have small couplings: 

£, A«1 , £ , A< 5 2 < 1 . (48) 

In fact this restriction arises naturally in the construc¬ 
tion of the effective theory [23]. As discussed follow¬ 
ing (8), the sigma model coupling g 2 = g 2 /T. Cor¬ 
rections at one loop order at one loop order are g 2 /T 
times this coupling; hence £/<? 2 and X/g 2 are rsj 9%I T , or 
£,A~(sf/T) 2 . 

Considered purely as a theory in two dimensions, it is 
possible to satisfy (48), and still have contributions at 
one loop order dominate over those at higher loop order. 
At one loop order, corrections are generically ~ g 2 , A, 
and £ times the quantity at tree level. At two loop order, 
corrections are then ~ g A , g 2 A, A 2 , and so on, times the 
quantity at tree level. Thus it is possible to take — for 
example — A and £ ~ g 3 . In this case, (48) is satisfied, 
since £ ~ A ~ g -C 1, but terms at one loop order, 
~ A and £, are still larger than two loop terms from the 
ordinary sigma model, ~ <? 4 . 


V. ONE LOOP EFFECTIVE ACTION 
A. What to Compute 

Having obtained the effective Lagrangian to quadratic 
order in (3, to one loop order the effective action involves 
the complete inverse propagator in the presence of the 
background field: 

S eff = i tr log (v ab + V? b di + . 


The trace involves summation over isospin indices and 
space-time momenta. 

In writing the inverse propagator, we take the liberty 
to freely integrate terms by parts. Terms in the inverse 
propagator can be categorized according to the number 
of times the derivative operator appears. Note that this 
corresponds to a derivative on the quantum held, (3, and 
not to a derivative on the background held, a. 

The term with two derivatives is a Laplacian in the 
background held, 

V ab = - d 2 S ab (50) 

+ 2 £ Qi sin 2 (|) di 6 ab 

+ 2 A (sin (^) a a ) di di (sin (^) a b ) . 

The hrst term on the right hand side, — <9 2 , is the usual 
Laplacian for a free held. 

The term free of the derivatives is a type of mass term, 

M ab = 2 A cos (|) <9 2 cos (|) S ab (51) 

+ f £o(a) (cos 2 (|) S ab -sin 2 (|) a°3 b ) . 

As for the second term in (49), we hrst consider a toy 
Lagrangian. Let a held /?“ interact through a derivative 
interaction with a background held V“ b : 

Cv = \(d i f3 a ) 2 + V‘/ b i3 a di/3 b . (52) 

We do not assume anything about the symmetry of V“ b 
in the isospin indices. The Laplacian in this background 
held is 

V£ b = - d 2 6 ab + i (vf di + di V ba ) . (53) 

The two terms in ~ Vi arise because the derivative can 
act either on the [3 which appears on the left, or that 
which appears on the right. For the last term, where the 
derivative acts to the left, it is permissible to integrate by 
parts. Doing so, the derivative acts either as an operator, 
or it acts on the background held: 

vf = ^d 2 s ab + ^ (v“ b - V ba ) di (d t v ba ) . (54) 

This result is general, as we have made no assumption 
about the order of the momenta. The result in (54) can 
also be derived more carefully, without any cavalier inte¬ 
gration by parts. One expands the original Lagrangian, 
(52), into terms which are symmetric and anti-symmetric 
in the isospin indices. This is done both for the back¬ 
ground held Vf b , and for /3 a dif3 b . After some algebra, 
one finds the same result as in (54). 

In the problem at hand, V“ b is given by 

K b = -i e “ bc (l + 2£ sin 2 (|)) A c 
+ £ sin(a) d a A b . 


(49) 


(55) 


Eq. (54) now contains one term that has one derivative 
on the quantum field, /3: 

Vt b = Vf~V^ a = -e obc (l + 2£ sin 2 (|))^ 

+ £ sin(a) ( ci a A b - a b A°;) , (56) 

and one term, — di V ba , which has no derivatives on /3: 

-di V-' 2 = - £ di ( sin (a) a b A“) • (57) 


We regularize integrals by dimensional continuation 
from two to 2 + e dimensions, but checked that the same 
results are obtained with Pauli-Villars regularization. We 
also derived all results using a mass to cutoff infrared di¬ 
vergences. We do not present the details of these checks, 
since they don’t affect the final results, but mention them 
to reassure the skeptical reader. 

B. Terms with no derivatives 


This is a type of mass term, and so we add it to that of The first example is trivial, the free energy of a mass 
(51), and define term: 


M ab = - £ $ ( sin(a) a b A?) 

+ 2 X cos 9 2 cos 5 ab (58) 

+ | ^o(<5) (cos 2 <5 Qb -sin 2 . 

The effective action is then 


- tr log (- d 2 + to 2 ) . 


(60) 


Computing in 2 + e dimensions, the ultraviolet divergent 
contribution to the effective Lagrangian is 


m 

~2 


d 2+e k 1 


(27r) 2+e k 2 


^ to . . 

72 ~ ~ 4^e ' (61) 


S eS = i tr log (V ab + V“ b d z +M ab ) , (59) 

where the term with one derivative now is like a coupling 
to a gauge field Vf b . 

Computing the complete effective action for this in¬ 
verse propagator would be most involved. However, our 
needs are much simpler. We do not want all terms in the 
effective action, but only the ultraviolet divergent terms 
which contribute to the renormalization of the original 
Lagrangian. For this, we only need to expand the effec¬ 
tive action to include terms with two derivatives of the 
background field, a. There is an infinite series of terms 
with higher numbers of derivatives of a, but these are 
ultraviolet finite. Terms with higher derivatives do have 
infrared divergences, but these are cutoff by masses gen¬ 
erated by the loop potential. 

The simplest term to compute is that from the mass 
term, A4 ab . Like the original Lagrangian, this is already 
of second order in derivatives of the background field, 
and so can be expanded directly. Terms involving the 
derivative operators, V“ b i di and V ab , are more subtle to 
compute. Some pieces of these operators are already of 
quadratic order in derivatives of a, but other pieces are 
not. Thus we have to take care, to ensure that we include 
all terms of quadratic order in the background field. 

Besides terms of quadratic order in derivatives of the 
background field, there are terms of zeroth order in 
derivatives. Such terms represent ultraviolet divergent 
renormalizations of the loop potential. In terms of the 
renormalization group, these contribute to anomalous di¬ 
mensions for the operators in the loop potential. We do 
not consider such terms in the present work, since the 
first thing to compute are the /3-functions. 

In the next three subsections, we discuss examples 
which illustrate how to compute each of these terms. We 
put all of the pieces together in Sec. (VI). 


This result is valid only to ~ 1/e, neglecting all contri¬ 
butions which are finite as e —> 0. This integral is needed 
to compute the contribution of A4 ab in (58). 


C. Terms with one derivative 


Since V“ b is linear in derivatives of the background field 
we have to expand it to quadratic order to obtain an ul¬ 
traviolet divergent term of quadratic order in derivatives. 
This contributes to the effective action, ^ tr logA ab , as 



V Qb 

3 


d 2+e k 

(2n) 2+e 


1 

167T6 



k l ki 
(fc 2 ) 2 ’ 


(62) 


The factors of i arise in going from coordinate to momen¬ 
tum space, <9i=> ik l . 


D. Terms with two derivatives 

We turn to the expansion of terms containing up to 
two derivatives. 

We start with the case where the background fields can 
be treated as constant. Then the Laplacian derived from 
(50) is just the function 

V ab = (l + 2/f sin 2 (I)) 6 ab 

+ 2 A sin 2 a a a b , (63) 

times the usual Laplacian for a massless field, — d 2 . 

For constant fields, it is easy computing the corre¬ 
sponding propagator. In fact, we can further simplify 
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our algebra by noting that in order to compute the (3- 
functions to one loop order, it suffices to use the form of 
this propagator which is valid only to linear order in £ 
and A. This propagator is 

A ab » ((l-2£ s in 2 (|))<5 ab 

— 2 A sin 2 a a a b ^ , (64) 

times the usual massless propagator, l/(— d 2 ). 

Notice that the restriction on the coupling constants 
in (48) is obvious in this form. If £ or A were not 1, 
they would alter the usual propagator of free field theory, 
and perturbation theory would be more complicated. 

As we explain in the next subsection, it is necessary to 
use this modified propagator in computing the effects of 
the terms with no and one derivatives of the background 
field. 

The problem is that in general, one is not allowed to 
treat the background fields as constant. To understand 
the problem, consider the toy model, 

C f = \{l + f)(d i p) 2 , (65) 

where / is some background field. We wish to compute 
the ultraviolet divergences for terms which are quadratic 
in derivatives of /. 

Besides such terms, there are also ultraviolet diver¬ 
gent terms which are constant in /; we ignore these in 
the present model, as renormalizations of the loop poten¬ 
tial. In two space-time dimensions, there are also terms, 
quadratic in derivatives of /, which arise through the con¬ 
formal anomaly [24] . Terms from the conformal anomaly, 
though, are ultraviolet finite, and so can be ignored in 
computing /3-functions. 

The Laplacian for this toy model is 

V/ = - d 2 +*d l fdi ■ (66) 

Terms of linear order in / can be neglected, since they 
only contribute to terms constant in /. The first terms 
which depend upon derivatives of / arise at quadratic or¬ 
der. We then expand the effective action, | tr log V/, to 
quadratic order. We go to momentum space, with p the 
momentum going through the background field, /. Then 
the term of quadratic order in the effective Lagrangian is 

1 fl )( f d2+ek ^ k ' +p )) 2 
4 J[P, \J (27r) 2 + e fc 2 (fc + p) 2 

The structure of the numerator is easy to understand. 
This is a one loop diagram, where the background field 
f(p ) couples to two quantum fields. The momenta at 
each vertex are k and k + p, with the coupling propor¬ 
tional to the product of these momenta. 

Since we only want the momentum dependent term, we 
can subtract off the value of the integral for p = 0. This 


) /(-*>)• ( 67 ) 


reduces a quadratically divergent integral to one which 
is merely logarithmically divergent, 

1 . / f d 2+e k (k ■ p) 2 - k 2 

~4 f[P) ^ J (27r) 2+e (fc 2 ) 2 


f{-p) . (68) 


Using the integral in (62), the contribution to the renor¬ 
malized Lagrangian is 


1 

167T6 


(5J) 2 , 


(69) 


where we have reverted to coordinate space. 

A slightly more involved toy Lagrangian, which most 
closely mimics our problem, is 


£f,h = + f) {m 2 + \idiih • /3)) 2 . (70) 


By similar computation, the ultraviolet divergent coun¬ 
terterms at one loop order are 


- T ^ e ((d i ff + 4(l-f)(d i h) 2 -2(d i h 2 ) 2 ) . (71) 


VI. /3-FUNCTIONS TO ONE LOOP ORDER 

With these results, we can compute the ultraviolet di¬ 
vergent contributions to the Lagrangian at one loop or¬ 
der. 

In doing so, it helps to recognize that this is all we 
want to do. We only want terms which renormalize the 
bare Lagrangian, and so only need keep terms ~ £o(a), 
£i(a), and £ 2 ( 0 1). At one loop order, we explicitly 
see that new interactions, as discussed at the end of 
Sec. Ill A, do appear. These include, for example, 
sin 4 (a/2)£o(a) ~ sin 2 (a/2)£i(a), and sin 2 (a/2)£ 2 (a). 
Such terms do have ultraviolet divergences, which con¬ 
tribute to the /3-functions for these couplings. We ignore 
these other couplings in the present work, considering 
the three couplings which we do include as representa¬ 
tive. Partial results about the renormalization of such 
couplings are given in Sec. VII, (84). 

Another important simplification is to recognize that 
terms at one loop order are one power of the coupling 
constant times the bare Lagrangian, (29). Remembering 
the definitions of couplings in (44), terms ~ g 2 times 
those in the original Lagrangian generate those ~ 1, £, 
and A; terms £ times the original Lagrangian generate 
contributions ~ £, £ 2 , and £ A; those A times the original 
Lagrangian generate terms ~ A, £ A, and A 2 . Thus at 
one loop order, we have to include all terms ~ 1, A, 
£ 2 , £ A, and A 2 . 

Of course, terms which are ultraviolet finite can be 
ignored completely. This includes all terms with more 
than two derivatives of the background field. 

The easiest contribution to compute is that from terms 
which are already of second order in derivatives, A4 ab in 
(58). Even so, we have to recognize that we cannot use 
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(61) directly, but must use the propagator in a back¬ 
ground field. Fortunately, as we already have terms with 
two derivatives, we can use the propagator for a constant, 
background field in (64): 

- — A ab M ab , (72) 

47T£ 


which equals 


1 

7re 



+ i(2£ + 3£A + 3£ 2 )£i (73) 

+ i (3 A +8^ + 16 A £ + 4 A 2 ) £ 2 j . 


Of course implicitly, all of the £’s are functions of the 
background field, a. 

The next contribution is from terms which are of first 
order in the derivatives, Vf fc in (56). Again, we cannot 
use the free propagator, but must use the propagator in a 
constant, background field. To one loop, this contributes 
the ultraviolet divergent terms as in (62), 


1 

167T6 


)c ycd da 


which is 


4(i£„ + i(A- 2 r) (£»-£,)) 


(74) 


(75) 


The final contribution is from terms with two deriva¬ 
tives. Using the toy model of (70), where we identify 


/ = 2 £ sin 2 (|) , h = 2 sin (|) a , (76) 

then from (71), we find that the ultraviolet divergent 
contributions to the one loop renormalized Lagrangian is 

-(-JaA) + 7 C A Ci (77) 

7 re \ 8 4 

+ I(a- 2 P + 4A 2 )a) • 


The sum of the ultraviolet divergent counterterms at 
one loop order, (73), (75), and (77), is 


Cct = 


+ 

+ 


\ (2£- A + 4CA + 5p)A 
i (3 A + £ 2 + 8£A + 4A 2 ) 


(78) 


We write the renormalized Lagrangian as 

r _r xMc | *aA 

Lren ~2Z g g2 L ° + Z*g* Ll + Z* g* 


(79) 


The renormalization constants for the three couplings are 
Z g , A, and Z\. In 2 + e dimensions, the couplings g 2 , £, 
and A all have dimensions of /A, where g is a renormal¬ 
ization mass scale. 

The renormalization constants are fixed by requiring 
that the ultraviolet divergences cancel in the sum of the 
counterterm Lagrangian, £ ct in (78), and the renormal¬ 
ized Lagrangian, C ren in (79) [24]. At one loop order, 
this determines the renormalization constants to be 


A 


1 

47T6 

1 

47T6 
1 

1 + “- 

47TC 


1 + 
1 + 


(. 9 2 - 3£ - A) 



1H —6a) , 
145-6 A -y) 


(80) 


The inverse powers of the coupling constant appear wor¬ 
risome; they certainly do not arise in the /3-functions of 
other nonlinear sigma models [24] . 

The /3-function for the coupling g 2 is given by 

P(g 2 ) = eff2 /( 1 _e ^ logZs ) ’ ( 81 ) 


and similarly for the other couplings. 

Using this definition, we find that the /3-functions have 
a standard form: 

Pig 2 ) = e g 2 + ^ ( _ 9 4 + 3 g 2 £ + g 2 A) + ... , (82) 

PiO = e £ + {— g 2 A + 11 f; 2 + 6 £ A) + ... , 

/3(A) = e A + (+ g 2 A + £ 2 + 14 £ A + 6 A 2 ) + ... . 


In the limit where £, A <C g 2 , the /3-functions reduce to 

Pig 2 ) ~ e g 2 - ^ g 4 + ■ ■ • > ( 83 ) 

Pit) ~ ^ 9 2 A + ■ ■ • , 

/3(A) ~ e A + — g 2 A + ... 

47T 

These are the principal results of our paper. We now 
discuss the properties of these /3-functions in two dimen¬ 
sions, where e = 0. 

As discussed in Sec. Ill A, the leading term in the 
/3-function for the coupling g 2 , Pig 2 ) ~ — g 4 /(47r), is 
the same as for the 0(4) nonlinear sigma model [24, 26] 
(accounting for the somewhat unconventional normaliza¬ 
tion of our coupling constant). As we compute in the 
limit where £ and A are <C g 2 , this term dominates those 
~ +<7 2 ? and ~ +g 2 A in (3(g ). As the leading term is of 
negative sign, the coupling g 2 is, inescapably, asymptot¬ 
ically free. 

The diligent reader might wonder why we didn’t save 
ourselves much effort, and directly compute the /3- 
functions in the relevant limit, where £ and A <tC g 2 . 
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In fact, this is what we did first. We computed the (3- 
functions for a field which is constant in isospin space, 
didi = 0, where it is much simpler to compute. As seen 
from (34), however, this only gives the /3-function for the 
sum of the two couplings, £ + A. After some labor, we 
found that this /3-function vanishes to the requisite order, 
~ g 2 £ and ~ g 2 A. From (83), we see that this is because 
in /3(£) and /3(A), the only terms of this order are ~ g 2 A. 
These terms are equal and of opposite sign, so that up 
to terms of higher order, they cancel identically. 

Thus in order to know the /3-function for £ + A, it is 
necessary to determine the nonleading terms in the /3- 
functions, as in (82). An important question is if it is 
possible to obtain an asymptotically free theory in all 
three couplings. Notably, since we take both £ and A to 
be g 2 , the theory is sensible for either sign of each 
coupling. When £ ~ A, however, the theory is clearly 
not asymptotically free: the dominant term in the /3- 
functions for both £ and A is ~ g 2 A. But as each term 
appears with opposite sign, (83), whatever sign we take 
for A, one coupling is asymptotically free, and the other, 
infrared free. 

Since A dominates the /3-functions, one can also con¬ 
sider the limit where we tune A to be much smaller than 
£, with A < ( <C j 2 . In this case, however, the full cal¬ 
culation, (82), show that the dominant term in /?(£) and 
/3(A) are each ~ +£ 2 . Each of these terms has a positive 
sign, and so produces infrared freedom, regardless of the 
sign of £. 

VII. CONCLUSIONS 

In this paper we computed the /3-functions for the sim¬ 
plest extension of a 0(4) nonlinear sigma model, which 
is a TV = 2 matrix model with three coupling constants. 
To leading order in weak coupling, the result is extremely 
simple, (83). In general, we find that there is no region 
of parameter space in which all three coupling constants 
are asymptotically free, (82). 

We conclude by discussing unpublished work by one 
of us [30]. Both of these calculations take a background 
field which lies along a single direction in the Lie algebra. 
As discussed in Sec. IV A, this ansatz is not adequate 
to give complete information about the /3-functions. As 
we shall see, however, these limited results indicate that 
when TV > 2, all matrix models include couplings which 
are not asymptotically free. 

As discussed in Sec. II, there is an infinity of couplings 
which we neglected. Consider a background field which 
lies along one direction in isospin space, = 0. To rep¬ 
resent the neglected couplings, we add to the Lagrangian 
the term 

(^+T Sin2 ” (f) {9ia)2 ' (84) 

For n = 1, this corresponds to (34), with the coupling 
Hi = £ + A. When n > 2, H„ is proportional to a sum of 


several, independent coupling constants. As in Sec. IIIB, 
in expanding about a trivial background, with a =>• g(3, 
this term first contributes to a 2(n + 1) point of /3’s. It is 
for this reason that we normalize this term by an overall 
factor of l/(< 7 2 ) n+1 , so that the 2(n + 1) point function 
has coupling H n , independent of g. 

The /3-function for H ra is found to be ~ +(n 2 — 1 )g 2 H n 
[30] . This is valid to ~ g 2 H n , and neglects subleading 
terms H 2 , etc. This vanishes for n = 1, in agreement 
with (83). When n > 1, though, the /3-function for H„ is 
infrared free. 

Calculations were also done for TV > 3. Following Sec. 
II, there are couplings, analogous to g 2 , £, and A, for all 
TV > 3. Going through the same calculations as for TV = 
2, it is straightforward to determine the leading terms in 
the /3-functions. The leading term in the /3-function for g 2 
is ~ —TV// 4 , so this coupling is, of course, asymptotically 
free [26]. When TV > 3, it is found that the /3-function 
for £ + A vanishes to leading order, ~ // 2 £ and ~ g 2 A 
[30]. This is precisely analogous to the results for TV = 2, 
and suggests that the form of (83) is valid for arbitrary 
TV, with /3(£) = —/3(A) ~ +g 2 A, up to corrections ~ £ 2 , 
~ £A, and ~ A 2 . That is, at least one of the two couplings 
is always infrared free. We also expect that when TV > 3, 
the infinite series of couplings to higher order, as in (84) 
for TV = 2, includes infrared free couplings. 

In summary, these partial calculations strongly sug¬ 
gest that the /3-functions for TV > 3 are very similar to 
those for TV = 2. Generically, while the dominant cou¬ 
pling g 2 is asymptotically free, there are always subdomi¬ 
nant couplings which are infrared free. We do not expect 
any region of parameter space in which all couplings are 
asymptotically free. 

We contrast this to the behavior of ordinary nonlinear 
sigma models [24, 25, 26, 28, 29]. These models only in¬ 
volve a single coupling constant, which is asymptotically 
free in two dimensions. In such theories, the symmetry is 
unbroken in two dimensions, and the theory is always in 
a massive phase. Above two dimensions, an ultraviolet 
stable fixed point appears, with g 2 ~ e. The appearance 
of this ultraviolet stable fixed point is assumed to be re¬ 
lated to a phase transition in which the symmetry breaks. 
The properties of the transition, as determined from this 
ultraviolet stable fixed point, working up in 2 + e dimen¬ 
sions, is believed to reflect the same universal properties 
as for a linear sigma model, working down from 4 — e 
dimensions [24]. 

There is a greater variety of phase transitions possible 
in matrix models, which is reflected in their Lagrangians. 
Matrix models have potentials, and we can couple them 
to gauge fields. In two dimensions, the terms in the po¬ 
tential, and the gauge coupling, all have dimensions of 
mass squared, and so dominate in the infrared limit. 

Conversely, since these terms have positive mass di¬ 
mension, they do not affect the ultraviolet behavior, in 
either two or 2 + e dimensions. From present results, and 
those of [30], it appears that the /3-functions are very 
similar for all TV > 2. 
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This is unlike the results of lattice simulations, which 
find that the order of the deconhning transition changes 
as a function of N. For theories in 2 + 1 dimensions, 
the transition appears to be of second order for N = 2, 
3, and possibly 4, and of first order for N > 5 [2], For 
theories in 3 + 1 dimensions, the deconhning transition 
is of second order for N = 2 [3], and of first order for 
N > 3 [4, 5, 6, 7]. 

This demonstrates that the connection between the ul¬ 
traviolet behavior of /3 -functions in 2 + e dimensions, and 
phase transitions in the infrared limit, is not as immedi¬ 
ate in matrix models, as it is in ordinary sigma models. 
Certainly, before we can understand what happens non- 
perturbatively in the infrared limit, we must first under¬ 
stand perturbation theory in the ultraviolet limit. This 
was the goal of the present work. We hope that it pro¬ 
vides an impetus to better understand the rich structure 
of phase transitions possible in matrix models. 
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IX. APPENDIX 

We use this opportunity to make a comment about 
the correlation functions of Polyakov loops at large N. 
Consider the two-point function of Polyakov loops, 

(tr L(a;)^tr L(0)) ~ exp(— m.i\x\) , x —> oo, (85) 

i 

where the summation is over all states which contribute 
to this correlation function. For a SU(N) gauge the¬ 
ory without quarks, the symmetry broken is Z{N ). As 
N —> oo, this becomes 17(1). By Goldstone’s theorem, 
when U{ 1) breaks, there must be an associated massless 
particle. At large but finite N, there is a light, almost 
Goldstone particle; as the mass of an almost Goldstone 
particle is proportional to the parameter which breaks 
the symmetry, which in the pure glue theory is 1/N 1 2 , 
then tUq ~ 1/-/V 2 , or mo ~ 1/iV. This almost Goldstone 
mode is separate from the usual, perturbative excitation 
associated with the Debye mass. 

This is trivial in three dimensions, but one might won¬ 
der if it still holds in two dimensions, since then a con¬ 
tinuous symmetry cannot break. We suggest, however, 
that at large N, the two-point correlations of (85) will 
look similar in either two and three dimensions. This is 
because at large N, all connected correlation functions 
are suppressed by powers of ~ 1/N 2 . Thus it is possible 
to have a massless Goldstone boson at infinite V, since 
then it is also non-interacting. 


[1] B. Svetitsky and L. G. Yaffe, Nucl. Phys. B 210, 423 
(1982), and references therein. 

[2] M. Gross and J. F. Wheater, Z. Phys. C 28, 471 (1985); 
J. Christensen, G. Thorleifsson, P. H. Damgaard and 

J. F. Wheater, Phys. Lett. B 276, 472 (1992); Nucl. Phys. 

B 374, 225 (1992); M. Teper, Phys. Lett. B 313, 417 
(1993); C. Korthals Altes, A. Michels, M. A. Stephanov 
and M. Teper, Phys. Rev. D 55, 1047 (1997), 

[arXiv:hep-lat/9606021] ; J. Engels, F. Karsch, E. Laer- 
mann, C. Legeland, M. Lutgemeier, B. Petersson and 
T. Scheideler, Nucl. Phys. Proc. Suppl. B 53, 423 (1997), 
[arXiv:hep-lat/9608099] ; P. Bialas, A. Morel, B. Peters¬ 
son, K. Petrov and T. Reisz, Nucl. Phys. B 581, 477 
(2000), [arXivdiep-lat/0003004] ; ibid. 603, 369 (2001), 
[arXiv:hep-lat/0012019] ; A. Hart, B. Lucini, Z. Schram 
and M. Teper, J. High Energy Phys. 040 (2000) 0006, 
[arXiv:hep-lat/0005010] ; G. V. Dunne, I. I. Kogan, 
A. Kovner and B. Tekin, J. High Energy Phys. 032 (2001) 
0101, [arXiv:hcp-th/0010201]; I. I. Kogan, A. Kovner 
and B. Tekin, J. High Energy Phys. 062 (2001) 0105, 
[arXiv:hcp-th/0104047]; B. Lucini and M. Teper, Phys. 
Rev. D 66, 097502 (2002), [arXiv:hep-lat/0206027] ; 

K. Petrov, [arXiv:hep-lat/0211028] ; P. de Forcrand 
and O. Jahn, Nucl. Phys. Proc. Suppl. B 129, 709 
(2004), [arXiv:hep-lat/0309153] ; K. Holland, M. Pepe 
and U. J. Wiese, Nucl. Phys. B 694, 35 (2004), 
[arXiv:hep-lat/0312022] ; P. Bialas, A. Morel and B. Pe¬ 


tersson, ibid. 704, 208 (2005), [arXiv:hep-lat/0403027] ; 
K. Holland, J. High Energy Phys. 023 (2006) 0601, 
[arXiv:hep-lat/0509041] ; J. Liddle and M. Teper, PoS 
LAT2005, 188 (2006), [arXiv:hep-lat/0509082] . 

[3] J. Engels, J. Fingberg, K. Redlich, H. Satz and M. Weber, 
Z. Phys. C 42, 341 (1989). 

[4] G. Boyd, J. Engels, F. Karsch, E. Laermann, C. Leg- 
eland, M. Lutgemeier and B. Petersson, Nucl. Phys. B 
469, 419 (1996), [arXiv:hep-lat/9602007] . 

[5] B. Lucini and M. Teper, J. High Energy Phys. 050 

(2001) 0106, [arXiv:hep-lat/0103027] ; B. Lucini, M. 
Teper, and U. Wenger, Phys. Lett. B 545, 197 (2002), 
[arXiv:hep-lat/0206029] ; J. High Energy Phys. 061 
(2004) 0401 [arXiv:hep-lat/0307017] ; B. Lucini, M. Te¬ 
per and U. Wenger, Nucl. Phys. B 715, 461 (2005), 
[arXiv:hep-lat/0401028]; H. Meyer and M. Teper, J. High 
Energy Phys. 031 (2004) 0412, [arXiv:hep-lat/0411039] ; 
M. Teper, [arXiv:hep-th/0412005] ; B. Lucini, M. 
Teper, and U. Wenger J. High Energy Phys. 
033 (2005) 0502, [arXiv:hep-lat/0502003] ; J. 

Kiskis, [arXiv:hep-lat/0507003] ; B. Bringoltz and 
M. Teper, Phys. Rev. D 73, 014517 (2006), 

[arXiv:hep-lat/0508021]. 

[6] B. Bringoltz and M. Teper, Phys. Lett. B 628, 113 
(2005), [arXiv:hep-lat/0506034] . 

[7] K. Holland and U. J. Wiese, [arXiv:hep-ph/0011193] ; 
K. Holland, M. Pepe and U. J. Wiese, 



13 


[arXiv:hep-lat/0309008]; Nucl. Phys. Proc. Suppl. 
B 129, 712 (2004), [arXiv:liep-lat/0309062] ; K. Holland, 
P. Minkowski, M. Pepe and U. J. Wiese, Nucl. Phys. B 
668, 207 (2003); [arXiv:hep-lat/0302023] ; Nucl. Phys. 
Proc. Suppl. B 119, 652 (2003), [arXiv:hep-lat/0209093] ; 
M. Pepe, PoS LAT2005, 017 (2006) [Nucl. Phys. Proc. 
Suppl. B 153, 207 (2006)], [arXiv:hep-lat/0510013] . 

[8] M. Laine, [arXiv:hep-ph/0301011] ; J. P. Blaizot, 
E. Iancu and A. Rebhan, [arXiv:hep-ph/0303185] ; 

U. Kraemmer and A. Rebhan, Rep. Prog. Phys. 67, 
351 (2004) [arXiv:hep-ph/0310337] ; J. O. Andersen 
and M. Strickland, Annals Phys. 317, 281 (2005), 
[arXi v: hep-ph/0404164]. 

[9] A. Peshier, B. Kampfer, O. P. Pavlenko and G. Soil, 

Phys. Rev. D 54, 2399 (1996); A. Peshier, B. Kampfer 
and G. Soff, Phys. Rev. C 61, 045203 (2000), 

[arXiv:hep-ph/9911474] ; Phys. Rev. D 66, 094003 (2002), 
[arXiv:hep-ph/0206229] ; [arXiv:hep-ph/0212179] ; 

A. Peshier, Phys. Rev. D 63, 105004 (2001), 

[arXiv:hep-ph/0011250]; A. Peshier, Phys. Rev. D 
70, 034016 (2004), [arXiv:hep-ph/0403225] ; A. Peshier 
and W. Cassing, Phys. Rev. Lett. 94, 172301 (2005), 
[arXivdiep-ph/0502138]. 

[10] E. V. Shuryak and I. Zahed, Phys. Rev. C 70, 

021901 (2004), [arXiv:hep-ph/0307267] ; Phys. Rev. 

D 70, 054507 (2004), [arXiv:hep-ph/0403127] ; 

[arXiv:hep-ph/0406100] ; G. E. Brown, C. H. Lee, 
M. Rho and E. Shuryak, Nucl. Phys. A 740, 171 
(2004), [arXiv:hep-ph/0312175] ; J. Liao and E. V. 
Shuryak, [arXiv:hep-ph/0508035] ; Phys. Rev. D 73, 
014509 (2006), [arXiv:hep-ph/0510110] ; E. V. Shuryak, 
[arXivdiep-ph/0510123]. 

[11] K. Kajantie, M. Laine, K. Rummukainen and 

M. E. Shaposhnikov, Nucl. Phys. B 458, 90 

(1996), [arXiv:hep-ph /9508379] ; ibid. 503, 357 

(1997), [arXivdiep-ph/9704416] ; K. Kajantie, 
M. Laine, J. Peisa, A. Rajantie, K. Rummukainen 
and M. E. Shaposhnikov, Phys. Rev. Lett. 79, 
3130 (1997), [arXiv:hep-ph/9708207] ; K. Kajantie, 

M. Laine, A. Rajantie, K. Rummukainen and 

M. Tsypin, .J. High Energy Phys. 011 (1998) 9811, 
[arXivdiep-lat/9811004] ; K. Kajantie, M. Laine, 
K. Rummukainen and Y. Schroder, Phys. Rev. Lett. 
86, 10 (2001), [arXiv:hep-ph/0007109] ; K. Kajantie, 
M. Laine, K. Rummukainen and Y. Schroder, Phys. Rev. 
D 67, 105008 (2003) [arXiv:hep-ph/0211321] ; J. High 
Energy Phys. 036 (2003) 0304, [arXiv:hep-ph/0304048] ; 
A. Hietanen, K. Kajantie, M. Laine, K. Rummukainen 
and Y. Schroder, J. High Energy Phys. 013 (2005) 0501, 
[arXiv:hep-lat/0412008] . 

[12] P. Giovannangeli, Phys. Lett. B 585, 144 (2004), 

[arXivdiep-ph/0312307] ; [arXivdiep-ph/0410346] ; Nucl. 
Phys. B 738, 23 (2006), [arXiv:hep-ph/0506318] . 

[13] M. Laine and Y. Schroder, J. High Energy Phys. 
067 (2005) 0503, [arXiv:hep-ph/0503061]; Proc. Sci. 
LAT2005 (2006) 180, [arXiv:hep-lat/0509104]. 

[14] D. J. Gross, R. D. Pisarski and L. G. Yaffe, Rev. Mod. 
Phys. 53, 43 (1981); N. Weiss, Phys. Rev. D 25, 2667 
(1982); K. Enqvist and K. Kajantie, Z. Phys. C 47, 
291 (1990); V. M. Belyaev and V. L. Eletsky, Z. Phys. 
C 45, 355 (1990) [Sov. J. Nucl. Phys. 51, 168 (1990)]; 

V. M. Belyaev, Phys. Lett. B 254, 153 (1991); A. 
Gocksch and R. D. Pisarski, Nucl. Phys. B 402, 657 
(1993). 


[15] T. Bhattacharya, A. Gocksch, C. Korthals Altes and 

R. D. Pisarski, Phys. Rev. Lett. 66, 998 (1991); 
Nucl. Phys. B 383, 497 (1992), [arXiv:hep-ph/9205231] ; 

S. Chapman, Phys. Rev. D 50, 5308 (1994), 

[arXiv:hep-ph/9407313] ; J. Wirstam, Phys. Rev. D 
65, 014020 (2002), [arXiv:hep-ph/0106141] ; P. Gio¬ 
vannangeli and C. P. Korthals Altes, Nucl. Phys. B 
608, 203 (2001), [arXiv:hep-ph/0102022]; ibid. 721, 1 
(2005), [arXiv:hep-ph/0212298] ; ibid. 721, 25 (2005), 
[arXiv:hep-ph/0412322] ; D. Diakonov and M. Oswald, 
Phys. Rev. D 68, 025012 (2003), [arXiv:hep-ph/0303129]; 
Phys. Rev. D 70, 016006 (2004), [arXiv:hep-ph/0312126] ; 
Phys. Rev. D 70, 105016 (2004), [arXiv:hcp-ph/0403108] ; 
E. Megias, E. Ruiz Arriola and L. L. Salcedo, Phys. 
Lett. B 563, 173 (2003), [arXiv:hep-th/0212237] ; Phys. 
Rev. D 69, 116003 (2004), [arXiv:hep-ph/0312133] ; 
[arXiv:hep-ph/0412308] ; J. High Energy Phys. 073 (2006) 
0601, [arXivdiep-ph/0505215]. 

[16] T. Banks and A. Ukawa, Nucl. Phys. B 225, 145 
(1983); R. D. Pisarski, Phys. Rev. D 62, 111501 (2000), 
[arXiv:hep-ph/0006205] . 

[17] A. Dumitru and R. D. Pisarski, Phys. Lett. B 504, 282 
(2001), [arXiv:hep-ph/0010083] ; ibid. 525, 95 (2002), 
[arXiv:hep-ph/0106176] ; Phys. Rev. D 66, 096003 
(2002), [arXiv:hep-ph/0204223] ; A. Dumitru, O. Scav- 
enius, and A. D. Jackson, Phys. Rev. Lett. 87, 
182302 (2001), [arXiv:hep-ph/0103219] . O. Scavenius, 
A. Dumitru and J. T. Lenaghan, Phys. Rev. C 66, 
034903 (2002), [arXiv:hep-ph/0201079] ; R. D. Pisarski, 
[arXiv:hep-ph/0203271] , and references therein; A. Du¬ 
mitru, D. Roder and J. Ruppert, Phys. Rev. D 70, 074001 
(2004), [arXivdiep-ph/0311119]. 

[18] P. N. Meisinger, T. R. Miller, and M. C. Ogilvie, Phys. 

Rev. D 65, 034009 (2002), [arXiv:hep-ph/0108009] ; 
P. N. Meisinger and M. C. Ogilvie, Phys. Rev. D 65, 
056013 (2002), [arXiv:hep-ph/0108026]; I. I. Kogan, 
A. Kovner and J. G. Milhano, J. High Energy Phys. 017 
(2002) 0212, [arXiv:hep-ph/0208053] ; K. Fukushima, 
ibid. 68, 045004 (2003), [arXiv:hep-ph/0303225] ; Phys. 
Lett. B 591, 277 (2004), [arXiv:hep-ph/0310121]; 

Y. Hatta and K. Fukushima, Phys. Rev. D 
69, 097502 (2004), [arXiv:hep-ph/0307068]; 

[arXiv:hep-ph/0311267]; P. N. Meisinger, M. C. Ogilvie 
and T. R. Miller, Phys. Lett. B 585, 149 (2004), 
[arXiv:hep-ph/0312272] ; A. Mocsy, F. Sannino and 
K. Tuominen, Phys. Rev. Lett. 91, 092004 (2003), 
[arXiv:hep-ph/0301229]; J. High Energy Phys. 044 
(2004) 0403, [arXiv:hep-ph/0306069]; Phys. Rev. 

Lett. 92, 182302 (2004), [arXiv:hep-ph/0308135] ; 

[arXiv:hep-ph/0401149] ; Jour. Phys. G 30, S1255 

(2004), [arXiv:hep-ph/0403160] ; F. Sannino and 

K. Tuominen, Phys. Rev. D 70, 034019 (2004), 

[arXiv:hop-ph/0403175] ; P. N. Meisinger and M. C. 
Ogilvie, [arXiv:hep-ph/0409136]; M. Schaden, Phys. 
Rev. D 71, 105012 (2005), [arXiv:hep-th/0410254], 

[19] O. Kaczmarek, F. Karsch, P. Petreczky, and F. Zantow, 

Phys. Lett. B 543, 41 (2002) [arXiv:hep-lat/0207002] ; 
S. Digal, S. Fortunato, and P. Petreczky, 
[arXiv:hep-lat/0211029] ; O. Kaczmarek, F. Karsch, 
P. Petreczky and F. Zantow, Nucl. Phys. Proc. 
Suppl. B 129, 560 (2004) [arXiv:hep-lat/0309121] ; 

O. Kaczmarek, S. Ejiri, F. Karsch, E. Laermann 
and F. Zantow, Prog. Theor. Phys. Suppl. 153, 
287 (2004), [arXiv:hep-lat/0312015] ; P. Petreczky 



14 


and K. Petrov, Phys. Rev. D 70, 054503 (2004), 
[arXiv:hep-lat /0405009]. 

[20] A. Dumitru, Y. Hatta, J. Lenaghan, K. Orginos and 
R. D. Pisarski, Phys. Rev. D 70, 034511 (2004) 

[arXiv:hep-th/0311223]. 

[21] O. Aharony, J. Marsano, S. Minwalla, K. Papadodimas, 
and M. Van Raamsdonk, Adv. Theor. Math. Phys. 8, 603 
(2004), [arXiv:hep-th/0310285] ; Phys. Rev. D 71, 125018 
(2005), [arXiv:hep-th/0502149]. 

[22] A. Dumitru, J. Lenaghan, and R. D. Pisarski, Phys. Rev. 
D 71, 074004 (2005), [arXiv:hep-ph/0410294], and refer¬ 
ences therein. 

[23] R. D. Pisarski, [arXiv:hep-ph/0608242] . 

[24] J. Zinn-Justin, “Quantum Field Theory and Critical Phe¬ 
nomena” (Clarendon Press, Oxford, 2002). 

[25] R. D. Pisarski, Phys. Rev. D 20, 3358 (1979). 

[26] A. McKane and M. Stone, Nucl. Phys. B 163, 169 (1980). 

[27] D. Friedan, Phys. Rev. Lett. 45, 1057 (1980); Annals 


Phys. 163, 318 (1985). 

[28] A. M. Polyakov and P. B. Wiegmann, Phys. Lett. 
B 131, 121 (1983); P. B. Wiegmann, Phys. Lett. B 
141, 217 (1984); P. Wiegmann, ibid. 142, 173 (1984); 
E. Ogievetsky, P. Wiegmann and N. Reshetikhin, Nucl. 
Phys. B 280, 45 (1987); J. Balog, S. Naik, F. Nie- 
dermayer and P. Weisz, Phys. Rev. Lett. 69, 873 
(1992); T. J. Hollowood, Phys. Lett. B 329, 450 (1994), 
[arXiv:hep-th/9402084] ; V. A. Fateev, V. A. Kazakov 
and P. B. Wiegmann, Nucl. Phys. B 424, 505 (1994), 
[arXiv:hcp-th/9403099] . 

[29] P. Fendley and H. Saleur, [arXiv:hep-th/9310058] ; 
P. Fendley, J. High Energy Phys. 050 (2001) 0105, 
[arXiv:hep-th/0101034] ; M. Caselle and U. Magnea, 
Phys. Rep. 394, 41 (2004), [arXiv:cond-mat/0304363] , 
and references therein. 

[30] M. Oswald, unpublished. 


